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ABSTRACT 


The in-line and transverse forces acting on cylinders of 
1.0 inch to 2.5 inches diameter placed near a plane wall ina 
harmonically oscillating flow have been measured. The drag 
and inertia coefficients C4 and Cra for the in-line force and 
the lift coefficients Cin and Cit for the transverse forces 
away and toward the wall, respectively, have been determined 
as a function of the relative gap e/D between the wall and the 
cylinder and the period parameter Vt/D- The relative gap 
ranged from 0.0143 to unity and the period parameter from zero 
to 35. In the subcritical range where these experiments have 
been performed, the effect of the Reynolds number, which ranged 
from 4,000 to 30,000, was coun to be secondary and certainly 
obscured by the excellent correlation provided by the period 
parameter. 

Tne frequency of the transverse force was found to range 
from one to five times the frequency of oscillation of the har- 


monic motion, but a dominant frequency of twice the oscilla- 


tion frequency was seen over most of the values of V/D- 
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NOMENCLATURE 


amplitude of the motion 

average drag coefficient 

maximum lift coefficient away from the boundary 
maximum lift coefficient toward the boundary 
average inertia coefficient 

diameter of test cylinder 


normal distance between the boundary and the cylinder 
surface 


instantaneous total force acting on the test cylinder 
drag force acting on the test cylinder 
lift force acting on the test cylinder 
water depth 

wave number (27/1) 

length of test cylinder 

wave length 

period of oscillation 

time 

horizontal fluid particle velocity 
instantaneous velocity 

maximum velocity 

vertical fluid particle velocity 
vertical position of the test cylinder 
percent error 

fluid kinematic viscosity 

fluid density 


wave frequency (217/T) 
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ZI. INTRODUCTION 


A. EXPERIMENTAL JUSTIFICATION 

With the current energy deficit becoming of vital impor- 
tance, offshore oil production assumed great priority, and 
pipe lines became the most feasible mode of liquid transporta- 
tion to the shore. This, in turn, necessitated the deter- 
Mination of the wave and/or current induced forces on submerged 
pipes placed at or near the ocean bottom. The problem, how- 
ever, 1s not unique to the pipe lines. 

Underwater cables for communication and Sensing purposes 
and all submerged structures, be they offshore oil platforms 
or oceanographic research instruments, are subjected to wave 
induced forces. 

The proximity of pipe to a solid boundary, such as the 
ocean bottom, further complicates the problem, and, unless 
the submerged pipe lines are anchored to a buoy system or 
buriea, the bottom effects must also be considered. Even 
buried pipe lines, subject to scouring can become uncovered 
and have their support dug out from under, leaving them sen- 


Sitive to wave induced forces. 


B. BACKGROUND THEORY 
The in-line force acting on a cylinder immersed in a time- 
invariant flow may be expressed, according to Morison and his 


co-workers [6!} as: 





_ . ———— 2 2 
Se 1/2LC Dev +C_ aD Via / at (1) 


in which Ca and oa represent time-invariant drag and inertia 
coefficients, respectively. The first term on the right-hand 
Side of this equation represents the velocity-squared depend- 
ent drag force and the second term the acceleration-dependent 
inertia force. 

A simple dimensional analysis shows that the time- 


averaged drag and inertia coefficients for a cylinder immersed 


in a harmonically oscillating uniform flow, represented by 
V=-V_cos 2ie7 < 


depend on the period parameter Viql/D and the Reynolds number 
V »p/Y- In addition to the in-line force, the cylinder is 
subjected to an alternating transverse force due to separa- 
tion and vortex shedding. The proximity of a solid boundary 
can alter significantly the magnitudes of both the in-line and 
transverse forces and requires the determination of the wall 
proximity effect on the drag, inertia, and lift coefficients. 
In fact, it is the determination of this effect that prompted 
the present study. 

A harmonically oscillating uniform flow represents only 
approximately the oscillatory motion induced by waves about 
a submerged cylinder. Nevertheless, the wave-induced forces 
may be predicted with sufficient accuracy through the use of 
drag, inertia, and the lift coefficients obtained with har- 


monic flow, provided that the velocities and accelerations 
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about the cylinder are either measured or calculated through 
the use of an appropriate wave theory. 

For small amplitude waves, Airy'‘s theory may be used to 
evaluate the velocities and accelerations of a given time and 
depth. The limitations of this theory are well-known and will 
not be represented here. 

The horizontal and vertical components of velocity may be 
WEitten as [1l], 

u=Ao(cosh k(h+z)/sinh kh) cos 6 (2) 

w=Ao (sinh k(h+z)/sinh kh) sin 6 (3) 
in which A represents the wave amplitude, k the wave number 


defined by 27/1, h the water depth, z the vertical position 


of the cylinder, o the wave frequency given by a and 
O=kx-ot. 
The components of acceleration are given by 
py 6t=-2ie - (cosh k(h+z)/sinh kh) sin 6 (4) 
Sew meme - (sinh k(h+z)/sinh kh) cos 6 (5) 


The in-line force may be expressed through the combination of 


equations (1), (2), and (4) to yield, 


F=A0°LC_ D2 cosh k(h+z) sin 86 
4 sinh kh 
2 2 
“AC, LoD (cosh _k(h+z)" | cosot | cosot (6) 
(sinh Bh) - 


The transverse force cannot, however, be calculated ina 


Similar manner through the use of equations (1), (3), and (5). 
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It is fundamentally the eddy behavior that determines the 
transverse force. Thus, the lift resulting from a nonlinear 
motion is not susceptible to estimation by superposition of 
individual responses to each harmonic acting in isolation. 

The lift coefficients in the present study are evaluated 
as follows: 


Cia (maximum transverse force away from the wall in a cycle) 
0.5 pDLV_* (7) 


Cio (maximum transverse force toward the wall ina Evcre) 
0.5 DLV_* (8) 


Force coefficients based on other time-invariant magnitudes 
such as the root-square values could have been easily obtained. 
It will suffice here to remark that a simple Fourier analysis 
for the in-line force and the peak value for the transverse 
force proved adequate for all purposes. 


Cc con- 


The determination of the coefficients CaCl Cine 1T 


stitute the essence of the present study for a cylinder placed 
within the proximity of a plane wall in a harmonically oscil- 


lating flow. 


C. RELATED RESEARCH 

Yamamoto, et al., in 1974, studied wave forces on cylinders 
near a plane boundary in the range of Reynolds numbers from 
2,000 to 30,000 and the surface period parameter from about 
0.3 to 3. The total force was considered as the sum of 


components due to water particle velocity squared (lift and 
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drag forces), and due to water particle acceleration (inertia 
force). They have correlated the force coefficients Ci and 


Cs with the relative distance of the cylinder from the bound- 
ary. 

The vertical and horizontal hydrodynamic forces were 
obtained for unseparated flow over the submerged cylinder by 
using the method of double images and applying Blasius' theorum. 

The results indicated that the proximity of a plane 
boundary modified both the lift and the inertia coefficients. 
However, these results are applicable only to unseparated 
flows where drag is negligible compared to inertial forces. 
Within a narrow range and relatively low values of the period 
parameter, it was shown that the maximum horizontal force oc- 
curred at zero crossings of the surface wave, indicating 
that wake dependent drag force is negligible and the horizon- 
tal force is composed mainly of the inertia force. The verti- 
cal force, for larger values of the relative distance was 
unsinusoidal with large downward forces (toward boundary) at 
the surface wave crests and small upward forces at surface 
troughs. 

In 1967, a theory for simple shear flow across a circular 
cylinder in proximity to a plane wall, was developed by Arie 
and Kiya [2]. Theoretically, they showed that a force acts 
downward (negative uplift) for uniform velocity, but a force 
acting upwards (away from the plane) exists for flow having a 


velocity gradient. Their experimental data failed to enforce 


1S 





this theory, but did show that lift on the cylinder increased 
as the clearance between the cylinder and wall decreased. A 
wind tunnel was used for their experimentation. 

Studies on the effects of ocean currents on pipes anchored 
just above the ocean floor were conducted by Wilson and Cald- 
well [12] in 1970, and the lift and drag coefficients were 
determined. The frequency of the forces on the pipe was 
found to be effected by the eddy-shedding frequencies. 

In 1971, a study by Grace [3] on the effect of the clear- 
ance of pipe above a flume and the effect of orientation of 
the pipe to the wave crest was made. With clearance ranges 
from 1/32 inch to 1-1/2 inches and wave periods of 2 to 6 
seconds in 3 foot deep water, he found that the horizontal 
force normal to the pipe was not effected by bottom clearance 
and that this force decreased as the pipe became perpendicular 
to the wave fronts. He also found that the "vertical force" 
decreased as the bottom clearance decreased. 

In much of the work conducted on wave forces, correlations 


between Care and the Reynolds number V p/ vs where v is 


qa! 
the kinematic viscosity, have been sought with little success. 
Wiegel, Beebe, and Moon (1957) [11] found no relation between 
these parameters in studies of ocean wave forces on cylindri- 
cal piles. Jen [4], in 1968, also failed to correlate Cha or 
Ca with Reynolds number in a laboratory study of a 6 inch 
diameter pile. 

In 1958, Keulegan and Carpenter [5], investigating the 


inertia and drag coefficients on circular cylinders, found 
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a correlation between these force coefficients and a dimen- 
sionless parameter, V/D- Working in a rectangular basin 
with standing waves, they covered Vt/D values between 2.7 and 
120 over a Reynolds number range trom 4000 to 29,300. The 
lowest value of Cn and the highest value of Ca were found to 
occur ata V_t/D Of abowte 15. 

Sarpkaya and Tuter [8], in 1974, with an oscillating U- 
shaped channel and various cylinders, reconfirmed Keulegan- 
Carpenter data except for a small Spread in the upper regions 
ot Vial/D» where the Sarpkaya-Tuter data may be more accurate 
and reliable. An important finding of the Sarpkaya-Tuter work 
was that the lift or transverse force was as larage as the in- 
line force. 


The close correlation between C_, and C., and the period 
1 


d 
parameter Vqi/D is evidenced by the close fitting curves re- 
lating these parameters. The Ch curve of Kuelegan-Carpenter 
forms a lower envelope for ocean data gathered by Wiegal, Beebe, 
and Moon, while the Cy curve forms an upper limit for the same 
Gatea (11). 

The work on wave induced forces on cylinders has attracted 
much attention, but results seem quite scattered and some- 
times contradictory. Also, the effect of a plane boundary in 
proximity to the cylinder presents a new variable. It is the 
purpose of this thesis to analyze the forces on a cylinder in 


the proximity of a plane wall and to correlate the various 


coefficients with the period parameter Vt/D, and the relative 


ris 





distance between the cylinder and the wall, throuch the use 


of harmonic motion. 
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II. METHOD OF ANALYSIS 


To analyze the time dependent force in unsteady periodic 
flow, a Fourier analysis will be made, beginning with the 
Morison's equation given by 


P=1/2 C,eAV|V|+ Coy av/dt (9) 


The enclosure of one of the components of v- term in "absolute 
value" brackets accounts for the change in sign of the resist- 
ing force as velocity changes direction in an oscillating 
flow. 


For a circular cylinder, the equation reduces to: 


F=C_m/4LD*pdv/at + Sa LDpv|v| (10) 
2 


Representing the velocity of the harmonic motion by 


V=-V_cosot (11) 
m 


and inserting into equation (10), one has 


2 > ge C_ sinot - C,|cosot|cosot (12) 


2 
pLV,, D/2 VT 
Multiplying both sides by sin t and integrating, one has, 


27 :; 7) 
y PE simot dicot) = me w=-C_sinotsinot d(ot) 
eV a D2 I m 


27 
-of C,|cosot|cosot|sinot d(ot) (13) 
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which gives 


27 ; = 2 
f F sanot ot) = 1 Dd C tt) (14) 
o 
pLV_D/2 ae 
and 
LN ati Fsinot d(ot) 
a ce” 8! Te ia 
D pLV D 


Multiplying both sides of the force equation by cosot, and 
integrating, results in 


27 


o F Gesot d(ot) = J — C_Sinot cosot d(ot) - 
"pL _“D/2 Vint 
(16) 
27 
oui C,|cosot cosot cosot d(ot) 
27 
eel Ey cosot d(ct) = -8/3 Ca (17) 
pLv D/2 
Finally, 
at Pacosgusd(ot) 
Ca = -~3/4 Q — (18) 


L 
le) ve D 


For the analysis of these force coefficients a computer pro- 
gram was used (See Appendix A). The determination of ven 1s 


made by considering the wave amplitude, A, and the period, T, 


resulting in 


Ma 27TA/T (19) 
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Thus the coefficients Ca and C; may be written as, 


2 eae 








Cs = 27 . iF sin - Peace. ) (20) 
D LAp 

Ca = -377 3 EFcos (425) ate 1} (21) 
8oDTLA 


As an alternate method of computing the force coefficients, 
a least-Squares method was also employed. The method of least 
squares consists of the minimization of the error between the 
measured and calculated forces. Letting F represent the in- 
stantaneous measured force and Pa the force calculated through 
the use of equation (10), and writing 

E*= (F-F_)* (22) 


and dE*/ac_=0 and dE*/dC .=0, one has 


27 
C.= = af pees congt d(ot) (23) 


: 37 oDLV 
m 
and 
27 
_ An = ioe 9 Sint d (ot) (24) 
D a LD 


It should be noted that the Fourier analysis and the method of 
least-squares yield identical Ca values and that the Cy values 
Gautier only slightly. 

The error between the measured and calculated forces, 


particularly in the neighborhood of the maximum forces, may 
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be further minimized by choosing the square of the measured 
force as the weighting factor in the least-square analysis. 


Thus writing 
E°=F* (F-F,)* (25) 


and dE*/ac ,=0 and dE*/dc_=0, one has 


C = 9 34S 2- (26) 
ome 21 3° 3 
and 
7? ie = fof 
1™~ LAD ae 663° 3 
The functions fy are given by 
27 D 4 271 3 
£,=f F cos ot d(ot) , £,= JS F |cosot| cosot d(ot) 
°o ° 
277 5 
£,=of FP sinot cosot| cosot|d(ot) (28) 
om > 3 27 3 
£y=of Pesimegot diet) , f= pire crnot ‘d(ot) 


Equations (27) and (28) may be shown to reduce to equations 
(23) and (24) by replacing FP? in equations (28) by one and 
Carrying out the necessary integrations in which F does not 
appear. 

Each wave cycle was divided into 36 time intervals of 


At=2.86/36 seconds, (T=2.86). The force for each interval was 


found ana incorporated into the computer program to calculate 


20 





the appropriate coefficients through the use of the three 


methods or analysis cited above. 
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Tit.) SEXP ERIMENTAL EQUIPMENT 


The basic oscillating flow system consisted of a U-shaped 
vertical water channel with a channel cross section of 18 by 
20 inches shown in Figure i. Oscillations were created by 
pneumatic pressure applied at a closed end cf the channel. 

The closed end was opened by means of a slider crank mechanism 
uncovering a large exit and releasing the air pressure. The 
resulting oscillation was a near perfect harmonic as evidenced 
by the elevation and acceleration traces shown in Figure 2. 

The maximum amplitude of the oscillation was 11 inches and 
the natural damping of the oscillations was in the order of 
1/8 inch per cycle. The water level at its minimum oscilla- 
tion height was twenty inches above the test cylinder. The 
cylinders were manufactured out of piexiglass tubes or rods at 
desired diameters and at lengths approximately 1/16 inch less 
than the channel width. Self-aligning bearings were imbedded 
at each end of the cylinders and caution was taken before each 
experiment to confirm the clearance of the cylinder end from 
the channel walls. 

The lateral and in-line (with respect to the fluid flow) 
force measuring devices consisted of various cantilever beams 
mounted to the ends of the test cylinder. Eight piezoresistive 
strain gauges were mounted on each cantilever beam and properly 
waterproofed. 

These fSBGe transducers were repeatedly calibrated by 


hanging known loads at the midsection of the cylinder in the 








Bitz J) Schematic 7 ada of a 


and cylinder orientation 
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horizontal and vertical directions. Some of the transducers 
were comprised of two pieces of thin cantilever beams cut 
orthogonal to each other and were able to measure both the 
in-line and lateral force simultaneously. Other transducers 
were capable of measuring only the in-line or lateral force. 
However, being rotatable “90°, this transducer could first 
measure the in-line force and with a 90° rotation, measure the 
lateral forces. In many cases the in-line and lateral forces 
were measured at each end of the cylinder and compared with 
each other. In no case did the force curves deviate from each 
other more than 5%, indicating a fairly uniform response along 
the cylinder with respect to the resultant forces. 

Throughout the investigation, the monitoring ct the char- 
acteristics of the oscillations in the U-channel was of prime 
importance. Most of the difficulties in past determinations 


OF Cail ane Gy resu@ted@irom tie diffire@lty of producing a 


Su Il 

purely harmonic motion or from having to determine oscilla- 

tion characteristics indirectly. The U-tube provides a per- 
fectly sinusoidal oscillation and the instantaneous displéece- 
ment and acceleration are continuously recorded. The instan- 
taneous elevation in one leg of the channel was determined 
through the use of a capacitance wire connected to an amplifier- 
recorder system. Such wires have been used in the past to 
measure wave heights in open channels. The response of the 


wire was found through calibrations, to be perfectly linear 


within the range of oscillations encountered. 
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The instantaneous acceleration was measured by means of a 
differential-pressure transducer connected to two pressure 
taps at the mid-section of one of the channel walls. The 


acceleration was then calculated from 
Ap= psa, (29) 


where A P is the differential pressure, op the fluid density, 
s the distance between pressure taps, and a. the instantane- 
ous acceleration. 

The effect of pressure drop due to the viscous forces was 
found to be negligible. 

The displacement and acceleration are, of course, in vonase 
and may be used independently to calculate the velocity and 
displacement of the fluid, although displacement was actuaily 
used. The smoothness of the variation of acceleration and 
elevation traces shows the success in obtaining a purely 


harmonic motion as shown in Figure 2. 
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Vee Dee CUS omOnNeOr RESULTS 


The drag coefficients are shown in Figure 4 as a function 
of the period parameter Viat/Ps which also equals 2nA/D, for 
various values of e/D. The actual data points are shown only 
for one value of e/D, partly to simplify the presentation of 
the data and partly to give some idea of the variation of Ca 
with respect to the three methods of analysis used in its 
evaluation. Firstly, it 1s apparent erat all.three methods 
yield nearly identical resuits and that the data exhibit very 
little scatter. Secondly, the drag coefficient can reach 
values as high as 3.75. For e/D larger than unity, Ca values 
approach those found for V yt/D=@. 

Evidently, it is not possible to explain the complex 
Variations of Ca with Viat/D Since it is largely determined by 
separation effects. 

The inertia coefficients are shown in Figure 5 in a manner 
Similar to that for Ca: The values of om GOrresponding to 
Vint /D=0 are obtained from the potential theory [10]. The ex- 
perimental values approach in all cases those predicted by 
the potential theory as V jt /D70- Once again, the inertia co- 
efficients approach those obtained for the limiting case of 
e/D=~ [8] for e/D values larger than unity. 

The lift coefficients C (force away from the wall) are 


1A 


shown in Figure 7 and the lift coefficients C (force towards 


dik 


the wall) are shown in Figure 6. Also shown in Figure 6 are 
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lift coefficients calculated through the use of the potential 
theory [10] for Vyl/D=0- Aside from the fact that the lift 
coefficients can become very large for cylinders near a wall, 
Figs. 6 and 7 also show that the predictions of the potential 
theory for separated flow may be grossly in error. Figure 7 
shows that for a small gap (e/D=0.0143), large lift forces 
act on the cylinder in a direction away from the wall and Cit 
is zero for all values of Vil /P larger than 15. According 
to the potential theory, a large net force exists toward the 
wall when a small gap exists between the cylinder and the 
wall [10]. Thus, the predictions of the potential theory can 
be used only for extremely small values of Vit/D fOr whieh 
separation does not occur. Evidently, the occurrence of separa- 
tion depends not only on e/D and V»r/D but also on the 
Reynolds number. No special effort has been made to determine 
the maximum values of Vt/D for a given set of e/D values 
_and Reynolds number below which the predictions of the poten- 
tial theory may be applied. 

The practical significance of the occurrence of large 
lift forces both away from and toward the wall in a given 
cycle of oscillation is self evident. Such alternating forces 
could give rise to large vertical oscillations. These oscil- 
lations can cause scour and fatigue and lead to larger in- 
line forces than those calculated from equation (12) through 


the use of the coefficients from Figs. 4 and 5. [In fact, it 
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is not difficult to imagine the occurrence of coupled in-line 
and transverse oscillations with disastrous consequences. 

The frequency fT of the alternating transverse force was 
evaluated in terms of the frequency f of the fluid oscilla- 
tion and it was found, as done previously [8], that ft. is 
not a constant fraction of f and that several frequencies may 
occur during a given cycle. This seems to be rather reason- 
able because of the fact that in any cycle the magnitude of 
the velocity varies from zero to Mis and that there can hardly 
be a single Strouhal frequency. In spite of the multiplicity 
of the frequencies of vortex shedding, the examination of the 
force records have shown that the vortices are shed with sig- 
nificant regularity over most of the Viat/D values at a domin- 
ant frequency of f =2f. Other frequencies of the alternating 
lift force were f£=f,3f,4f, ame Site 

Two additional matters need to be discussed, namely, the 
error between the measured and calculated in-line forces and 
the importance cf the Reynolds number on the variation and 
correlation of the various force coefficients. 

The relative error in each cycle for each Viat/D as well 
as the error corresponding only to the maximum in-line force 
was calculated for each method of data analysis. The detailed 
results which will not be reported here for the sake of brev- 


ity have shown that } defined by 
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—— ee YE, (calculated) 


is within *15% and reaches its maximum in the neighborhood 

of Vzt/D=10 for the in-line force. This result may be taken 
as an indication of the fact that Morison's equation may be 
used to predict the in-line force acting on cylinders in the 
Vecimm@rty sat a wales in a hasmoniucally oscillating: fluid through 
the use of the Fourier averaged coefficients. 

As to the Reynolds number cCGefined by V/V s its role in 
wave force coefficients is not well understood. For Reynolds 
numbers below the critical value (a value which is expected to 
be considerably smaller than that corresponding to the steady 
flow about a cylinder, partly due to unsteady nature of the 
flow and partly due to the presence of the vortices both up- 
stream and downstream of the cylinder), the Reynolds number 
appears to play some role in the variation of Cy in the 
vicinity of the intermediate values of VjT/D- For V At/D about 
10, only one or two vortices form on the downstream side of 
the cylinder during a given cycle. As previously observed by 
Sarpkaya [9], when the fluid accelerates rapidly at high 
Reynolds numbers, vorticity is slow to diffuse and eneesrore 
accumulated rapidly in the growing separation bubble behind 
the cylinder. Although this bubble of trapped vorticity 
soon reaches unstable proportions, the growth of the bubble 


and, hence, of the resulting vortices, are so rapid that the 
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vortices beccme much larger than their quasi-steady-state 

size. Furthermore, the relatively slow diffusion of vorticity 
(largely due to turbulence rather than molecular diffusion) 
enables the vortices to retain a larger fraction of the 
vorticity shed from the boundary layers as the vortices move 

a small distance downstream before the flow is reversed. Dur- 
ing this period, the drag on the cylinder rises above its 
steady-state value, both because the surface area covered by 

the separation bubble is temporarily greater than it would 
ultimately become in steady flow and because the enlarged 

wake distorts the outer flow, reducing the base pressure below 
its steady value. Thus, the effect of the Reynoids number, 
which is nearly absent in steady flow within the range of 
Reynolds numbers (4,000 to 30,000) discussed herein, is primari- 
ly due to the differences in the rate of shedding and diffu- 
Sion of vorticity between the steady and periodic flow. How- 
ever, the turbulent nature of diffusion and the excellent corre- 


mation of C, with Vial/D nearly obscure the effect of the 


d 
Reynolds number, at least below the critical Reynolds number 
range. Near or above the critical Reynolds number range, how- 
ever, the effect of the Reynolds number may become more pro- 
nounced due to different reasons. It must be kept in mind 

felra t Ca and Ch are time averages and the error between the 


measured and calculated forces through the use of these time- 


invariant coefficients is not uniform and depends on V t/D- 
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Similar arguments hold true for the transverse forces. The de- 
pendence of Cin and Cit on Vial/D is much stronger and tends 


to obscure the weaker dependence on the Reynolds number. 


by 





z . V. CONCLUSIONS 


Although the in-line and transverse forces acting on 
cylinders near a wall were measured for sinusoidally oscillat- 
ing flow rather than for wave-type fluid motion, the results 
presented herein enable a number of conclusions concerning 
the latter type of flow to be drawn. 

Morison's equation may be used to predict the in-line 
force acting on such cylinders or deeply submerged pipes (no 
free surface effects) subjected to local wave motion through 
the use of the Fourier averaged coefficients. The drag and 
inertia coefficients approach those predicted by the poten- 
tail theory as V{t/D approaches zero. 

The transverse forces, like the in-line forces, are 
strongly influenced by flow separation. Not only the magni- 
tude but also the direction of the forces differ significant- 
ly from those predicted by the potential theory. For intermed- 
1ate values of the period parameter, potential theory under- 
estimates the transverse forces whereas for large values of 
Viai/Ds it grossly overestimates them. 

Additional work with regard to the practical pipeline 
problems is needed to ascertain the effects of free surface, 
nonlinear waves, internal currents, roughness, and other 
environmental conditions. The data presented herein could 


Porm the basis for Stch extensions. 
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APPENDIX A 
(IN-LINE FORCE COEFFICIENT PROGRAM) 
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